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ABSTRACT 



Upcoming surveys with the Herschel Space Observatory will yield far-IR photometry of large samples of young stellar objects, which 
will require careful interpretation. We investigate the color and luminosity diagnostics based on Herschel broad-band filters to identify 
and discern the properties of low-mass protostars. We compute a grid of 2,016 protostars in various physical congurations, present 
the expected flux densities and flux density ratios for this grid of protostars, and compare Herschel observations of three protostars to 
the model results. These provide useful constraints on the range of colors and fluxes of protostar in the Herschel filters. We find that 
Herschel data alone is likely a useful diagnostic of the envelope properties of young stars. 
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1. Introduction 

A comprehensive theory of star formation is necessary for un- 
derstanding processes spanning the entire range of cosmic evo- 
lution, from the formation and evolution of galaxies to the 
formation of planets. Thus, observing the formation of stars 
is one of th e prime science obj ectives of the Herschel Space 
Observatory Pilb rat et al. (2010)| With its large (3.5 m) primary 
mirror, Herschel is able to survey a substantial fraction of the 
current star formation in the Galaxy with unprecedented effi- 
ciency and sensitivity. The result will be an explosion of far-IR 
and sub-mm data on star formation regions. Herschel is partic- 
ularly suited for carrying out these surveys because the spectral 
energy distributions (SEDs) of protostars peaks in the far-IR and 
sub-mm region. A large number of guaranteed- and open- time 
key (legacy) Herschel programs are, thus, focused on star for- 
mation regions. The se include PACS (pho todetector array cam- 
era & spectrometer, 



& imaging reciever, 



Poglitsch etal. 2010| ) and SPIRE (spectral 
Griffin et al. 2010b Imaging surveys of the 



Gould belt iAndre et al. 201 OK and the Herschel Orion protostar 
survey (HOPS) and DIGIT programs, among others. 

These surveys will detect hundreds of thousands of proto- 
stars. For our purposes, protostars are low-mass stars in the ear- 
liest stages of their evolution. They are less than a few x 10 5 yrs 
old and are still surrounded by an infalling envelope. The SEDs 
of these protostars will be complex mixture of their intrinsic 



* Herschel is an ES A space observatory with science instruments pro- 
vided by Principal Investigator consortia. It is open for proposals for 
observing time from the worldwide astronomical community 
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emission coupled with reprocessing and emission from their im- 
mediate environment. Radiative transfer models are an impor- 
tant tool in understanding and interpretting the observed flux 
densities of protostars. Model SEDs can provide crucial con- 
straints on the most important physical parameters of a proto- 
stellar system. Because protostars can encompass a wide rage 
of parameter space it is often necessary to compare observa- 
tions against a large number of different models to derive a rea- 
sonable solution and to understand the degeneracies involved. 
IRobitaille et al. 20061 constructed a grid of 20,000 individual 
models for SED fitting; however, these models used an ism 
model for the envelope, which is not accurate at far-IR wave- 
lengths. This is because reprocessing by dust is especially impor- 
tant at HerscheVs wavelengths as this reprocessed light forms the 
bulk of the protostars' SED. In this contribution, we have gen- 
erated a grid of 20,160 SEDs from 2,016 individual protostaiQ 
models focused primarily on the parameter space appropriate for 
young stars, and for the Herschel broad-band filters. 



2. The model grid and flux predictions 

We use the Monte Ca rlo radiative transfer code developed by 
|Whitney et al. (2003)} which is adapted to include the dust mix- 
ture properties from Tobin et al. (2008) to calculate the model 
SEDs. For convenience, we briefly summarize the main com- 
ponents here. The code comprises four physical components: 
(i) a central object, (ii) A flared disk, (iii) A rotating collaps- 
ing envelope, and (iv) a bipolar cavity. These components are 



There are 10 different inclination angles for each of 2,016 protostars 
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illustrated in Fig.[T]. The star+disk system emits energy packets 
("photons") that are absorbed or scattered by the disk and/or en- 
velope, subsequently determining the emergent SED of the sys- 
te m. Our adopted dust mo del uses the methodology described 
by |D'Alessio et al. (200 1)| with addi tional optical constants for 
graphite from |Draine & Lee (1984)| The dust grains used here 
are a mixture of graphite, silicate and water ice and their number 
density has a power-law distribution with grain size to the power 
-3.5 and include large (1 micron radius) grains; the implied gas- 
to-dust ratio is 133. 

The code was used to produce a grid of 20,160 model SEDs 
using 2xl0 7 photons for the Monte Carlo computations, and 
by using various configurations for the main components listed 
above. We restricted the variable parameter set to: the infall rate 
(M gnv ), the maximum envelope size (R en v,max), the cavity open- 
ing angle (6c), the cavity shape exponent (b), the accretion rate 
(Mdtsk), the size of the disk (Rdisk,max), and the inclination an- 
gle of the system to the observer's line of sight (/). These are the 
critical parameters that predominantly determine the flux density 
prediction for Herschel PACS and SPIRE bands (56-610 //m). 
The remaining para meters have fixed v alues that are typical for 
low-mass protostars (ITobin et al. 20081) . 

One important characteristic of the code is that the disk ac- 
cretion rate is used in the model grid to vary the luminosity of the 
protostellar systems; the central star is fixed at 1 L©. This is rea- 
sonable since we do not know the relative contribution of pho- 
tospheric and accretion luminosities in protostars. However, the 
spectrum of the central object's emission does not largely affect 
the emergent spectrum of the protostar in the far-IR since nearly 
all emission is reprocessed. Thus, dust opacities are a more crit- 
ical and relevant factor for our analysis. 

Table [T] summarizes the relevant model parameters and 
their values in o ur grid. We list only select parameters (see 
ITobin et al. 20081 for a complete list). The properties of the cen- 
tral object are perhaps the most crucial of the fixed parameters 
in that changing the parameters for central object may lead to 
degeneracies in the observed colors and fluxes. We discuss this 
further in Sect. 13.51 Other relevant limitations of the model are 
that dust opacities, which are known to vary across regions, are 
assumed to be constant, and that only the water ice opacities for 
the 3.6 yum feature are included in addition to dust. 




Cavity 



Envelope 



Fig. 1. Illustration of some of the major parameters used to 
model the protostar SEDs. 



In protostellar systems, model grids, such as the one 
discussed here, are most accurate and useful when multi- 
wavelength data are available to provide as many constraints as 
possible (see e.g. IFischer et al. 20101) . Nonetheless, PACS and 
SPIRE sample the peak of the protostar' s SED and we expect 
Herschel data alone to provide a relatively good estimate of the 
envelope properties in particular. 

We further investigate the effect of each of the variable pa- 
rameters on the observed fluxes. To illustrate how a particular 
model parameter affects the flux densities, we first selected a typ- 
ical protostar system, here defined to have the following values: 
M disk =l0- 7 M Q yr" 1 , R disk , max =l00 AU, M^-10" 5 M yr" 1 , 



Re 



,=10,000 AU, C =25 degrees, b = 1.5, and i = 63.3 de- 
grees. Then, we show how individual parameters affect the re- 
sults as they vary over the range listed in Table [T] while hold- 
ing all other parameters constant. The results are color-coded (as 
listed in the legend boxes in the upper left hand part of Figs.[2l-[6]) 
and shown as thick dashed lines. We show and discuss below re- 
sults from a representative set of correlations investigated from 
the full grid of model parameter range. We omit discussion on 
M disk , as this parameter and its use in the code is already dis- 
cussed above (see Sect. [2]) . 



2.1. Synthetic Photometry 

The code computes model SEDs for a given spatial extent (aper- 
ture) of the system. We use aperture radius of 20,000 AU. 
This collects the total integrated flux from the envelopes which 
are all smaller than the aperture. We convolved the resulting 
SEDs with the filter + detector transmission curves of PACS and 
SPIRE (kindly provided by their respective instrument teams). 
We present the results as flux density values at the central wave- 
lengths of the filters. The model flux es are calculated at th e dis- 
tance of Orion, adopted to be 420 pc (IMenten et al. 20071) . 

3. Results and discussion 

Figures [2]-[6] show the results for the 20,160 model SEDs calcu- 
lated for our grid for various combinations of flux densities and 
flux density ratios. The shaded dots shows the full range of pre- 
dicted flux densities for the entire set of 2,016 protostar systems 
as seen with 10 different inclination angles each. The parameter 
values used here sample a broad range of physical conditions in 
the protostars and, thus, provide a reasonable comparative range 
of expected flux densities and flux density ratios for protostars. 



3.1. Envelope Infall rate 

The mass infall rate strongly affects the far-IR flux and not sur- 
prisingly, the predicted flux density and their ratios typically 
vary by factors ~2 as a function of M env . This is straight forward 
as increasing the amount of infall simply increases the amount 
of cold material and thus the far-IR fluxes. We expect to be able 
to constrain the envelope infall rates with relatively good preci- 
sion with Herschel data alone, although, some knowledge about 
the distance to the source (or its luminosity) is required, or must 
be assumed. 



3.2. Envelope size 

The size of the envelope directly affects the amount of the cold- 
est dust present in the protostar. Thus, we note significantly more 
variation in the longer wavelength SPIRE filters (see Figs.[5]-[6]) 
than in the PACS filters. Increasing the outer radius of the en- 
velope preferentially creates a comparatively redder protostellar 
system. However, the affect of the R en v,max is still appreciably 
lower than M env . That's because the latter more directly con- 
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Parameter Description 

Central object & observer: 



Value(s) 



M,(M ) 
R, (Ro) 


Mass of the central object. 
Radius of the central object. 


0.5 
2.09 


r.(K) 

;' (degrees) 


Surface temperature of the central object. 
Inclination angle 


4000 
18.2,31.8,41.4,49.5,56.7,63 


Envelope Properties: 






M^CxlO^Moyr" 1 ) 


Mass infall rate 


1,5,7.5, 10,25,50,100 


Kenv,min \K*J 

R c (AU) 


Envelope outer radius 

Envelope inner radius 
Centrifugal radius 


5000, 10000, 15000 

42.75 

Equal to R d i S k,max 


Disk properties: 






M disk (Mq) 


Disk mass 


0.05 


M^^xlO^Moyr" 1 ) 

Rdiskjnax (AU) 


Disk accretion rate/Luminosity fl 
Disk outer radius 


1, 10, 100, 1000 
5, 100, 500 


K-disk,min 

a 

p 


Disk inner radius 

disk radial density exponent 

disk scale height exponent 


Set by dust destruction radius. 

2.125 

1.125 


Cavity properties 






6 C (degrees) 
b 


Cavity opening angle 
Cavity shape exponent 


5, 15, 25, 45 
1.5,2.0 



Notes. 

Parameters in bold are model grid variables. 

(a) See Sect.[2]for the coupling between M disk and luminosity 
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Fig. 2. Predicted results for the PACS 70 yum filter vs. the ra- 
tio of 70 to 160 yum flux density. The fluxes are reported for 
a protostar at the distance of Orion. The shaded region shows 
the full range in the 20,160 SEDs in the grid. The color-coded 
(inset) thick dashed lines illustrate the magnitude and direction 
of the parameter when all other parameters are held constant. 
We do this by choosing typical values (see text) and plotting the 
variation from only one parameter for that protostar. The dia- 
mond symbols show the locations of the 3 HOPS protostars (see 
Sect.r 



trols the envelope density and total amount of emitting material 
present in the envelope. 

3.3. Cavity opening angle 

The PACS fluxes appears to be most strongly affected by this 
parameter (see Figs. [2] and |4]). This is likely because the width 
of the cavity preferentially carves away the part of the envelope 



Fig. 3. Same as Fig.[2]but this time for the 3 parameters that ap- 
pear to produce little or no effect on the observed flux densities. 
This Fig. is a magnified portion of the range of values shown in 
Fig.[2j The effect of these parameters will not be visible if shown 
over the full range displayed in Fig. [2 



responsible for the 70 /mi flux. By contrast, there is almost no 
variation as a function of 6c in the SPIRE filters. 



3.4. Inclination angle 

The optical and near-IR emission from protostellar systems is 
dominated by the source itself and the scattered light from 
the circumstellar envelope. Thus, face-on models are typically 
brighter and bluer compared to those in which the central star 
is either wholly or partially obscured by the disk. At optical 
and near-IR wavelengths, these objects are quite often extremely 
faint or undetected even in unbiased photometry surveys due to 
high obscurataion. Figures [3] shows that in the PACS filters this 
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3.5. The central object 

We have concentrated on a single central star; however, the mod- 
els are applicable to all low mass stars. The observed fluxes and 
SEDs simply scale with the mass of the central object as deter- 
mined by the envelope density profile and the accretion luminos- 
ity equations. When the mass of the central star is changed, the 
envelope density remains the same but the infall rate changes. 
And, the mass accretion rate for a given luminosity scales in- 
versely with the central object mass. Thus, considerable overlap 
does exist in the observed fluxes for different mass protostars. 
This result further demonstrates that PACS and SPIRE are best 
suited for constraining the envelope and disk properties and pro- 
vide limited constraints on the central object itself. 



Fig. 4. The ratio of the PACS flux densities, 70 /im/100 //m vs. 
100yum/160yum. The symbols and lines are as shown in Fig. [2] 



3.6. The HOPS protostars sample 
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Fig. 5. The SPIRE 250 fim flux density vs. the ratio 
250 yum/350 yum. The symbols and lines are as in Fig.[2j 
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Fig. 6. The PACS (70yum/160yum) ratio vs. the SPIRE 

(250yum/350yum) ratio. 



effect becomes largely unimportant as the main emission is from 
the envelope. We find similar results for the other combinations 
of flux and flux density diagnostics (not shown). 



The HOPS program dFischer et al. 20101) imaged 4 protostars in 
Orion L 1641 region for which PACS 70 and 160 yum photom- 
etry is available for 3 sources. These stars are shown with the 
diamond symbol in Fig. [2] Based on this diagnostic alone, we 
can conclude that all 3 HOPS sources are much more luminous 
than our typical protostar of 1 L©. Only 3 sources with two-band 
photometry are in the HOPS field; nonetheless, the sources are 
not clustered but spread over a broad range of color and flux den- 
sity values in Fig. |2l implying that a broad range of values must 
also exist for the other parameters. These conclusions are indeed 
confirmed by detailed model SEP fits ove r the 1-870 yum range 
for the HOPS sources dFischer et al. 20101) . 

4. Conclusions 

The model grid presented here provides the expected flux den- 
sities and flux density ratios for 20,160 low-mass protostellar 
SEDs, useful for determining the phase space occupied by pro- 
tostar in the Herschel filters. Disentangling the contributions 
of the various physical components that make up a protostel- 
lar system is observationally challenging and requires multi- 
wavelength data to properly constrain the properties of the sys- 
tem. Nonetheless, preliminary results show that the fluxes in the 
Herschel PACS and SPIRE filters are dominated by reprocessed 
light and emission from the protostellar envelope. Herschel data 
alone is, therefore, expected to be able to discern envelope prop- 
erties with some precision provided that a reasonable estimate 
for the object's distance or its luminosity is available. The source 
of the emission (combination of central star and disk accretion) 
is relatively unimportant at Herschel far-IR and sub-mm wave- 
lengths. Thus, Herschel data is unlikely to constrain the mass of 
the central object, but sampling the peak of the protostars SED 
allows tight constrains on its total luminosity. 
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